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We investigate transport in several translationally invariant spin-1/2 chains in the limit ol high 
temperatures. We concretely consider spin transport in the anisotropic Heisenberg chain, the pure 
Heisenberg chain within an alternating field, and energy transport in an Ising chain which is exposed 
to a tilted field. Our approach is essentially based on a connection between the evolution of the 
variance of an inhomogeneous non-equilibrium density and the current auto-correlation function at 
finite times. Although this relationship is not restricted to the case of diffusive transport, it allows 
to extract a quantitative value for the diffusion constant in that case. By means of numerically 
exact diagonalization we indeed observe diffusive behavior in the considered spin chains for a range 
of model parameters and confirm the diffusion coefficients which were obtained for these systems 
from non-equilibrium bath scenarios. 
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I. INTRODUCTION 



Although transport in low-dimensional quantum 
systems has intensively been investigated theoretically 
in the past years, there still is an ongoing interest 
in understanding the transport phenomena in such 
systems, including their temperature and length scale 

^^pj,^j^j^^^l. 2. 3.4. 5. 6. 7. 8. 9 . 10. 11. 12. 13. 14. 15. 16. 17.18. 19. 20. 21. 22. 23 . 

Those works have often addressed a qualitative classifi- 
cation of the occurring transport types into ballistic or 
normal diffusive behavior and, in particular cases, the 
crucial mechanisms which are responsible for the emer- 
gence of diffusion have been studied. In this context the 
role of non-integrability and quantum chaos is frequently 
discussed as an at least necessary condition^ i^°i^^i^^i^° . 
Significant theoretical attention has been devoted to spin- 
1/2 chainc^-^-^-^-^-^-^-^-^ -^°--^^-^^'^^-^^'-^^-^^'^^-^^-^^-^°-^^-^^-^'^ 
e.g. to the prominent anisotropic Heisenberg chain 
(XXZ model)i'^'^^'^>^>2'^>aia^i2JL3a4i5^i7i8^_ Most 

controversial appears the question whether or not the 
(finite temperature) transport in the pure Heisenberg 
chain is ballistiol^^^^^^i^^^^^. 



Even though there certainly is a large variety of differ- 
ent methods for the investigation of transport in quan- 
tum systems, we concentrate here on two of the main 
approaches in detail which are also most relevant in the 
context of the present work. The first approach may be 
classified as a direct one, since it is rather close to an 
experimental measurement setup: Within the theory of 
open quantum systemaS? the model of interest is coupled 
locally, e.g. at both ends of a spin chain, to reservoirs 
of the transported quantity, e.g. at different tempera- 
ture or chemical potentia^J^ii^i^Si^i^^i^i^i^^. Due to 
the coupling to reservoirs, a stationary non-equilibrium 
state eventually results for which all relevant expectation 
values such as the current and the spatial density profile 
of the transport quantity can be evaluated. Here, a van- 
ishing profile corresponds to ballistic behavior, whereas 



normal diffusive dynamics is associated with a strictly lin- 
ear profile, according to Fourier's law, respectively Pick's 
law, see, e.g. Ref. [ij. In the latter case a finite con- 
ductivity is simply given by the ratio of the current and 
the spatial density gradient of the transported quantity. 
This conductivity can be also understood in terms of a 
diffusion coefficient, since transport is driven by a gradi- 
24.2f"'^{) l^^ifeS |j]ji^^ considered model and not by an external 
lOTCc" *. One may therefore directly compare with the dif- 
fusion constant of a corresponding closed scenario, where 
transport is not induced by the coupling to baths but by 
an initially inhomogeneous non-equilibrium density^iiM. 

In order to simulate the influence of the baths the 
Liouville-von Neumann equation for the coherent evo- 
lution of the density matrix is routinely extended by an 
incoherent damping term, see, e.g. Ref. |33j. The deriva- 
tion of such a dynamical equation from a microscopic 
bath model is highly non-trivial and involves a combi- 
nation of various subtle approximation schemes, e.g. im- 
proper approximations may eventually lead to a mathe- 
matically correct but physically irrelevant quantum mas- 
ter equation (QME)^'^. One often intends to derive a 
proper QME of the Lindblad form^^ , because its special 
structure allows to apply the numerically efficient method 
of stochastic unraveling--?'-??^ e.g. spin-1/2 chains with 
about 16 sites become numerically tractable'^- . However, 
a recently suggested matrix product operator Ansatz has 
been shown to significantly increase this number of avail- 
able sites up to several dozensi^, i.e. finite size effects are 
drastically reduced and the extracted conductivity may 
be interpreted as a pure bulk property of the model, if 
it does not depend crucially on the concrete form and 
strength of the bath coupling, of course^^. 

Another approach for the investigation of transport 
in quantum systems is the Green-Kubo formula (KF) 
which was originally derived for electrical conductance 
by the use of linear response theor y^^'^^ . In that case 
the (frequency-dependent) electrical field is an external 
force which perturbs the system and the resulting current 
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of charge through the system is the response to this exter- 
nal perturbation. The KF as such gives the (frequency- 
dependent) conductivity as the hnear response coefficient 
in terms of a current auto-correlation function. The 
same approach is also used in the context of gradient- 
driven transport phenomena, e.g. for transport of energy 
or heat. In that case the current is driven by a much more 
complicated mechanis m^^i^^i^^i'^^i^^i'^^i'^"^ , especially since 
it can not be treated as a perturbation to the system. 
Nevertheless, simply by the replacement of the electri- 
cal current, e.g. by the energy current^, the KF is used 
for gradient-driven transport, too. But a rigorous jus- 
tification of this replacement remains still a conceptual 
problem, see Ref. [2^ or the overview paper |4l| (and the 
comprehensive literature which is cited therein). 

However, the KF is nowadays a standard 
method for the investigation of transport in spin 
chainsii^iii^i^iii^ii^iSSi^, not least due to its direct com- 
putability, once a finite piece of the considered system 
has been exactly diagonalized. As far as numerics is 
concerned, exact diagonalization is restricted to spin-1/2 
chains with at most 24 sites. But another difficulty 
arises, if the KF is evaluated on the basis of a finite piece: 
One often distinguishes between ballistic and normal 
diffusive transport by the notion of the Drude weight, 
essentially the conductivity's singular contribution at 
zero- frequency, see the review 0, for example. Whenever 
it is finite, the long-time behavior is expected to be 
ballistic, and whenever it vanishes, the dc- conductivity 
as the zero-frequency limit of the conductivity's regular 
non-singular part determines the long-time behavior. 
Since for any finite system the conductivity exclusively 
consists of delta peaks at different frequencies, there ap- 
parently are only singular contributions and it therefore 
is a difficult question how to extract and extrapolate the 
dc-conductivity from a finite system^. Moreover, the 
quantitative comparison of the resulting dc-conductivity 
(due to a possibly hypothetic external force) with a 
diffusion constant (from non-equilibrium bath scenarios) 
is hardly possible without the detailed knowledge about 
an Einstein relation between both quantities. 

In the present paper we do not intend to further dis- 
cuss the above mentioned conceptual and methodologi- 
cal problems which may come along with the KF. In- 
stead we will present a different but in a sense related 
approach which particularly is not concerned with the 
most of those problems. To this end we will firstly intro- 
duce a "typical" inhomogeneous non-equilibrium density 
in Sec. Ulland then connect the evolution of the variance 
of this density to a current auto-correlation function at 
finite times"^^. Remarkably, in the special limit of in- 
finitely long times this connection will be shown to yield 
a generalized Einstein relation which relates the diffusion 
constant to the dc-conductivity, i.e. as evaluated by the 
KF. Moreover, we will demonstrate that the great ad- 
vantage of the connection is given by its direct applica- 
bility at finite times and for finite systems as well. In this 
context we will suggest another concept for the analysis 



of data which is available from current auto-correlations 
functions. 

By the use of the suggested concept we will investigate 
in Sees. IllljfVl transport in several translationally invari- 
ant spin-1/2 chains in the special limit of high tempera- 
tures. We will concretely consider spin transport in the 
anisotropic Heisenberg chain, the pure Heisenberg chain 
within an alternating field, and energy transport in an 
Ising chain which is exposed to a tilted field. By means 
of numerically exact diagonalization we indeed observe 
strong indications for diffusive behavior in the considered 
spin chains for a range of model parameters and, what is 
more, we are able to quantitatively confirm the diffusion 
constants which were found for these systems from non- 
equilibrium bath scenarios in Refs. [Ulliilllll. Finally, 
we will close in Sec. lVII with a summary and a conclusion. 

II. CONNECTION BETWEEN VARIANCE AND 
CURRENT AUTO-CORRELATION FUNCTION 

In this Sec. we are going to introduce our approach 
to density-driven transport in translationally invariant 
quantum systems. To this end the next Sec. Ill Al firstly 
presents the pertinent definitions and exclusively de- 
scribes the general theory which eventually yields the 
basic Eq. ([13]) for the time-dependent diffusion constant 
!?(<). In the following Sec. lIIBl wc then motivate to eval- 
uate P(i) for finite times t and particularly illustrate the 
concept which will be used for the concrete spin chains in 
the subsequent Sees. Thus, the reader which is not pri- 
marily interested in the theoretical details may directly 
continue with Sec. IIIBI 

A. Diffusion Constant 

In the present paper we will investigate translationally 
invariant, one-dimensional quantum spin systems which 
are described by a respective Hamiltonian H. In those 
quantum systems we will consider an overall conserved 
transport quantity X, i.e. [_ff,Ar] = 0. This transport 
quantity and the Hamiltonian as well are decomposable 
into N formally identical addends i^, respectively 
corresponding to different positions, i.e. 

N N 

x = Y.^^, H = Y.K- (1) 

Thus, X ^ IS cL local density of the transported quantity 
X. Note that the may be defined on the positions of 
the /i^, in between, or both. The above decomposition is 
further done in such a way that Heisenberg's equation of 
motion for the local densities reads 

^ = « [H, x^] = I [h^,- , x^J + I [h^,+ , , (2) 

=j(M-i) =-j^L 
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where /i^- and represent those local addends of the 
Hamiltonian H which are located directly on the l.h.s., 
respectively r.h.s. of x^. This apparently implies a kind 
of locality. However, such a description can always be 
at least approximately enforced, if only interactions are 
reasonably short-ranged. For all quantum systems in the 
following Sees. IIIlllVl the description will be even exact, 
because interactions between nearest-neighbors are taken 
into account solely. As routinely done, the comparison 
with a continuity equation suggests the definition of a 
local current according to the scheme in Eq. see 
Ref. [25I, for example. This definition is consistent, if 



M+1) 



0, 



(3) 



where the latter holds, if X is globally preserved. The 
total current J is given by 



N 



(4) 



Once the above decomposition has been established, 
we can define a certain class of initial states p(0) which 
corresponds to some inhomogeneous, non-equilibrium 
density. To those ends let 



(5) 



denote the deviation of the local densities from their 
equilibrium average (a:^) — Trji^poq}, where pcq is any 
stationary equilibrium state, i.e. [H^ p^q] = 0. Then the 
initial state p(0) reads 



N 



p(0) = poq + X! Pcq 



(6) 



with some realization for the numbers Sf^. The factor 
is concretely given by 



(7) 



and therefore quantifies the equilibrium fluctuations of 
the transported quantity X. 

For the special initial state p(0) we now consider the 
actual expectation values 



It follows that 



d^{t)=Ti:{d^{t)p{Q)}. 



N N 

fi—i fi—i 



(8) 



(9) 



i.e. the sum 5 of the numbers 5^ in Eq. ^ determines 
the sum of the actual expectation values (t) . 



Of particular interest is the spatial variance W'^{t) of 
the d^{t). It is given by 



W^{t) 



S 

ti=i 



N 

E 



-I 2 



(10) 



If the dynamics of the (i^(t) was indeed generated by a 
discrete diffusion equation of the form 



dt 



d^{t)^V{t)[d^^i{t)~2d^{t)+d^+i{t)] , (11) 



then the evolution of this variance would read 
d 



dt 



W^t) = 2V{t), 



(12) 



as long as the d^{t) vanish at the ends of a chain (open 
boundary conditions) or are reasonably concentrated at 
a sector of a ring (closed boundary conditions). Even 
though Eq. pT|) implies Eq. p^ . the inverse direction 
generally is not true, of course. 

However, in Ref. |23 a connection of the form has 
recently been found directly from Eq. (jlOp . i.e. by firstly 
applying Heisenberg's equation of motion to Eq. (jlOp and 
subsequently manipulating the resulting equations. The 
time-dependent diffusion constant 2?(i) reads 



dt' C{t') 



(13) 



and is essentially given in terms of a time-integral over 
the current auto-correlation function 



C{t)=Tl{J{t)piqJ piq}, 



(14) 



i.e. C{t) = {J{t) J) in the limit of high temperatures 
(T — > 00). This limit will be considered throughout this 
work. 

Strictly speaking, the above diffusion constant T>{t) 
is restricted to the special initial state p(0) in Eq. 
of course. Nevertheless, due to its concrete form, this 
initial state represents an ensemble average w.r.t. to 
typicality^i^i^ or, more precisely, the dynamical typ- 
icality of quantum expectation values^. Thus, the over- 
whelming majority of all possible initial states with the 
same c?^(0) as p(0) is also expected to approximately 
yield the d^{t) corresponding to p(0), if the dimension 
of the relevant Hilbert space is sufficiently large, see es- 
pecially Ref. H^. The latter largeness is certainly ful- 
filled for all practical purposes. Or, in other words, 
the concrete curves for c?^(t) and thus for W'^{t) only 
slightly "fluctuate" around the curve of the ensemble 
average. In fact, for the single-particle quantum sys- 
tem in Refs. ,26,29,49 it has been demonstrated that 
Eqs. 1121)- (HH) correctly describe the dynamics for all pure 
initial states IV'(O)) which are not created explicitly in or- 
der to violate these equations, see Ref. [25ll49l . 

Eventually, let us relate the non-perturbative result 
for the diffusion constant V{t) to the standard result of 
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linear response theory for the dc-conductivity o'dd^'^. 
In the hmit of high temperatures this conductivity may 
be written as 

adc = |: ^°°dt'C(i'), (15) 

where (3 is the inverse temperature. The comparison of 
Eq. Hi]) with Eq. yields 



V{t oo) 



gdc 



(16) 



i.e. it leads to an Einstein relation. It is well-known that 
the integral in Eq. will diverge, whenever C(w), the 
Fourier transform of C{t), has a finite contribution at 
ljJ — 0, see Ref. 0, for example. As routinely done, we 
may hence define a Drude weight D, e.g. 



t~.oo t Ne^ 



B. Finite Time and Length Scales 



(17) 



If the total current is strictly preserved, i.e. [H, J] = 0, 
the diffusion constant I?(t) is completely governed by the 
Drude weight D, i.e. 'D{t) = D t iov all t. Consequently, 
transport is purely ballistic at each time, respectively 
length scale. Even if the total current does not represent 
a strictly conserved quantity, a non-zero Drude weight 
directly implies that the diffusion constant 'D{t) is still 
approximately given by the straight line Dt in the limit 
t ^ oo. However, at some finite time or length scale the 
diffusion constant X'(t) may nevertheless appear to be 
almost constant, i.e. Vit) « I?, as expected for diffusive 
behavior. Such a behavior requires that the Drude weight 
D is not relevantly large (on the corresponding scale), the 
finiteness of D by itself is not crucial in this context. In 
fact, it is well-known that the classification of transport 
types for a given quantum system generally is a concept 
which crucially depen ds on the considered time or length 
scale, see Refs. I25ll26l l29l for example. 

The above line of reasoning becomes also relevant in 
those situations where some model parameter A induces 
a transition from a finite towards a zero Drude weight 13, 
e.g. where D vanishes above some critical value of A, say 
Ac- Then at this critical value Ac a sharp transition from 
ballistic towards non-ballistic transport is to be expected 
at the infinite time scale. However, such a transition does 
not necessarily appear at some finite time, respectively 
length scale. Even if there was a transition, a sharp one 
would require at the critical value Ac a sudden jump of 
the Drude weight D from a finite and relevantly large 
number (on the corresponding scale) to zero. 

We hence do not concentrate merely on Drude weights 
and their finiteness, although those are also discussed, of 
course. Instead of that we focus on the diffusion constant 
'D{t), as defined in Eq. at finite times and for finite 
systems as well. According to Eq. (IT51) . we can evaluate 



the underlying current auto-correlation function directly 
in the time domain, i.e. C(i), instead of the frequency 
domain, i.e. C{lo). This way we are not concerned with 
the problems which may arise due to the fact that C{lij) 
is a highly non-smooth function for a finite system. We 
particularly exploit that the £'(i)-curve does not change 
at sufficiently short time scales any more, when the size 
of a system becomes large enough. Thus, at those short 
time scales, interesting signatures of the infinitely large 
system may already be extractable for a system with an 
accessible size. In the following Sees. IllUfVl this concept 
will be demonstrated in full detail. 



III. ISING CHAIN WITHIN A TILTED FIELD 
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FIG. 1: The diffusion constant ©(i), as given by Eq. (fT3|l . 
for energy transport in the Ising chain within a tilted field in 
the high temperature limit (T = oo). Insets zoom in ©(t) at 
short t. Parameters: J = —8, Bx = 6.75, = 4 (top) as 
well as Bz = Q (bottom). All curves are evaluated numerically 
(exact diagonalization) for chain lengths A'' = 8, 10, . . ., 16 
(arrows). The circles (top inset) represent additional data 
for N = 24, extracted from Ref. [2l| and computed by the 
use of approximative numerical integrators. The cross (top 
inset) indicates the conductivity ©bath = 11 from Ref. Rfa^. 
as obtained from a non-equilibrium bath scenario for the same 
set of parameters. 

In the present Section we will study a first example of 
a translationally invariant spin system. This system is 
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FIG. 2: The long-time increase D of the diffusion constant 
©(t), the Drude weight, for energy transport in the Ising 
chain within a tilted field at high temperatures (T = oo). 
Parameters: J = —8 and Bx = 6.75. All open symbols are 
evaluated numerically (by the use of exact diagonalization) 
for chain lengths A'^ < 16 and for z-components < 6. The 
filled square is extracted from Ref . [2l| and represents an upper 
bound for iV = 24. 



an Ising chain which is exposed to a, say, tilted magnetic 
field B. Concretely, its Hamiltonian reads {H = I)i^i2ii2^ 



N 



[ + ^ + ^ i^l+^l+i) > (18) 

where Bz and denote the z-component, respectively 
x-component of the total vector B = {Bx,QtBz). Here, 
one might think of some magnetic field which originally 
was in line with the z-direction and has been rotated 
about the y-axis with the angle a — aicia,n{Bx/ B^). 

In the Hamiltonian (jlSp the operators (i — x, y, z) 
are the standard Pauli matrices (corresponding to site 
/J,); J denotes the coupling strength; and N represents 
the total number of sites. Throughout this work we use 
periodic boundary conditions, i.e. + 1 = 1. 

Due to the presence of the tilted magnetic field, the 
only trivial symmetries are translational invariance and 
mirror symmetry^^ such that the whole Hilbert space is 
decomposable into 2 A^ decoupled subspaces with similar 
dimensions. For the present model the component of the 
total spin S = 1/2 Y^^=i ^i-i B-direction (and any other 
direction) is non-preserved. Therefore magnetization (or 
spin) is not a suitable transport quantity here. 

However, this model is appropriate to investigate the 
transport of energy, i.e. it allows to study the density 
dynamics of the local quantities /i^, see Eq. p^ . Note 
that these quantities contain contributions from both the 
Ising interaction and the Zeeman energy. The respective 



N 



Bx J 



and the factor e^, as defined in Eq. ([7]), reads 



16^"^ 



2Bi + 2Bi 



(19) 



(20) 



If the component B^ is identical to zero {a — 90°), the 
energy current (|19p is strictly preserved. Furthermore, 
the Hamiltonian is well-known to be integrable. But 
it can become quantum chaotic for Bz ^ 0, e.g. for the 
special set of parameters J = —8, Bz = 4, Bx ~ 6.75 
{a « 5Qo-j i3, 21^22 ^ Recently, for exactly this parameter 
set, a strong evidence for diffusive behavior has been 
found from non-equilibrium bath scenario o^'^i^^'^^ . 

Fig. [2 (top) shows the diffusion constant 2?(i), as given 
by Eq. p^ . for the above set of parameters. We observe 
that I?(t) increases within the correlation time of the 
underlying current auto-correlation function but, already 
on this time scale, essentially the Drude weight governs 
the overall shape of the curve, at least for small lengths 
N ^ 8. But when A^ is increased, the Drude weight D 
rapidly becomes smaller, i.e D decreases faster than a 
power-law, see Fig. [5] (squares). In particular there is no 
need to assume a finite and relevant value for D in the 
thermodynamic limit N oo. Already for A^ « 16 the 
Drude weight D is such small that the tendency of the 
diffusion constant I?(t) to gradually develop towards a 
horizontal line becomes visible. Nevertheless, only from 
those lengths which are available from numerically exact 
diagonalization a definite conclusion may still be vague, 
particularly conclusions about the time after which T>{t) 
possibly remains constant and about the constant value 
V, too. 

Fortunately, additional data is available in literature: 
In Ref. |2lj the current auto-correlation function has been 
evaluated also by the use of an approximative numer- 
ical integrator. Based on this data the time evolution 
of the diffusion constant I?(t) can be extracted^ for the 
length A^ = 24, see Fig. [T] (top inset). And in fact, for 
times which are larger than i « 10 we observe that V(t) 
takes on a constant value V « 10.5. The latter value 
further is in excellent agreement with the conductivities 
2?bath ~ 10.3 from Ref.|23and Pbath = H from Ref.lH^, 
as obtained therein for the same set of parameters from 
non-equilibrium bath scenarios. (In these works Pbath is 
denoted by k.) 

When the z-component is decreased from w 4 down 
to 0, diffusive behavior eventually breaks down towards 
ballistic transport, because the current becomes strictly 
preserved for Bz = 0- Unfortunately, it is hardly possible 
to give a critical value for this transition, simply due to 



the limited system sizes which are accessible by the use 
of exact diagonalization. However, such a value may be 
very unsharp for a continuous transition, i.e. if the Drude 
weight changes smoothly and does not jump suddenly 
from zero (or a finite, irrelevantly small number) to a 
finite, relevantly large number, cf. Sec. |TT1 

Contrary, when the z-component is increased from 
Bz ~ 4, the indications of diffusive behavior become 
much more pronounced for the accessible system sizes, 
cf. Fig. [1] (bottom). Respective diffusion constants may 
therefore be suggested for this parameter regime. 




FIG. 3: The diffusion constant 'D{t), as given by Eq. (|13p . for 
magnetization transport (spin transport) in the anisotropic 
Heisenberg chain (XXZ model) in the high temperature limit 
(T = oo). Insets zoom in ©(t) at short t. The curves are 
evaluated numerically (by the use of exact diagonalization) 
for the anisotropy parameter A = 0.5 and for chain lengths 
= 8 and 20. 



IV. ANISOTROPIC HEISENBERG CHAIN 
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In this Section we will investigate transport in the 
anisotropic Heisenberg chain (or XXZ model) as another 
and certainly more interesting example of a translation- 
ally invariant spin system. It is described by a Hamilto- 
nian of the form {h ^ Vp^L^^ 

N 

with the anisotropy parameter A. Independent from the 
concrete choice of A, the Hamiltonian (PT|) is integrable 
in terms of the Bethe Ansatz, see Ref. [jj, for example. 

In the presence of an external (uniform) magnetic field 
B one may add to the Hamiltonian ()2ip a Zeeman term 
of the form Hb = B S^. However, because the following 
investigation (at infinite temperature) will not depend on 
the concrete choice of _B, we set = for simplicity. 



FIG. 4: The diffusion constant V{t), as given by Eq. ITS]), for 
magnetization transport (spin transport) in the anisotropic 
Heisenberg chain (XXZ model) in the high temperature limit 
(T — oo). Insets zoom in ©(t) at short t. The curves are eval- 
uated numerically (by the use of exact diagonalization) for 
chain lengths = 8, 10, . . ., 20 (arrows) and for anisotropy 
parameters A — 1.0 (top), 1.5 (middle), as well as 2.0 (bot- 
tom). One cross (middle inset) indicates the conductivity 
©bath/J = 0.58 in Ref. [TJa, as found therein from a non- 
equilibrium bath scenario. Another cross (bottom inset) rep- 
resents the value a^d {P^ J) ~ 0.56 with the dc-conductivity 
Odc/iPJ) = 0.14 according to ReflH. 



As well-known, the Hamiltonian (j2ip is invariant under 
rotations about the z-axis, i.e. it commutates with . As 
a consequence H can be diagonalized within decoupled 
subspaces with dimensions N over M + N/2, where M 
is the quantum number w.r.t. S^. Moreover, due to the 
translational invariance as well as the mirror symmetry 
of (HH), the problem can be reduced further by a factor 
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FIG. 5: The long-time increase D of the diffusion constant 
T>{t), the Drude weight, for magnetization transport (spin 
transport) in the anisotropic Heisenberg chain (XXZ model) 
at high temperatures (T = oo). The symbols are evaluated 
numerically (by the use of exact diagonalization) for chain 
lengths A'^ < 20 and for anisotropy parameters A < 8. The 
line indicates a fit for the special case of A = 1.5. 

2N. 

Since represents a strictly conserved quantity, this 
model allows to investigate the density dynamics of the 
local quantities ct^/2, i.e. magnetization transport (or 
spin transport). The respective magnetization current 
is given by^iiii^ii^ 

N 



"f, = ' ^2^1^ K] = jK^l+i ~ ^l^l+i) (22) 

and the factor e^, as defined in ([7]), takes on the fixed 
value 1/4. 

Similarly, an energy current can be also introduced via 
the relation — i [ft,^, /i^+i]. But, since its commutator 
with the Hamiltonian (|2ip vanishes exactly, this current 
is strictly conserved such that energy transport is purely 
ballistic at each time- and length-scale (and for all finite 
temperatures, see Ref. Hol). 



Except for the special case of A = (XY model) the 
magnetization current is non-preserved. On that account 
it could be in principle possible to find diffusive behavior 
for non-zero anisotropics. Nevertheless, in the literature 
there is strong evidence that such a behavior is restricted 
to the regime A > 1 onljii^iii^i^ili^iiiii^iiiii^, most 
promising appear those A which are close to LSi^ii^ii^. 
Those expectations are also due to results on the Drude 
weight IJiiii^, e.g. Bethe Ansatz approaches suggest that 
D is finite for the regime A < 1 and zero for the regime 
A > 1^'-. Still controversial is the special case of A = 1, 
where D may be already zeroi or noli^. 

For completeness, we also show numerical results for 
the Drude weight in Fig. [SJ^ although respective data can 
be found already in Ref. 4, also obtained by the use of 
exact diagonalization and additionally extrapolated to 
^ oo. This data can be transferred directly to the 
present investigation, if the concrete values for the Drude 
weight are multiplied by a factor l/(e^7r) = A/it. Further 
results on Drude weights from exact diagonalization can 
be found in, e.g. Refs. [mH. 

A. Anisotropies A < 1 

Fig. [3] shows the diffusion constant X'(t), as defined in 
Eq. (|13p. for the anisotropy parameter A = 0.5. For TV = 
8 the Drude weight D determines the overall shape of the 
I?(f)-curve, i.e. on all time scales V^t) is very close to the 
straight line D t. There is only an insignificant deviation 
from this line at short times scales below the correlation 
time of the underlying current auto-correlation function, 
see Fig. [3] (inset). For N — 20, i.e. when the length of 
the chain is more than doubled, the X'(t)-curve remains 
practically the same. On that account it appears to be 
justified to assume a similar curve in the thermodynamic 
limit N oo. This assumption is consistent with the 
extrapolation of the Drude weight in Ref. d, see Fig. [5] 
(squares), too. We may therefore suggest purely baUistic 
behavior for the case A = 0.5. 

Note that Sirker et al. have recently presented results 
in Ref. [l3 which point towards a coexistence of diffusive 
and ballistic dynamics, where the ballistic distribution 
is small due to a small Drude weight. Even though we 
reproduce that the Drude weight is large for < 20 and 
does not depend significantly on TV, we can not exclude 
the possibility that the Drude weight eventually becomes 
small for N ^ 20, of course. 

B. Anisotropies A > 1 

In Fig. 3] (middle) we display the diffusion constant 
I?(i) for the anisotropy A = 1.5. We observe that 'D{t) 
increases at short time scales, i.e. within the correlation 
time of the underlying current auto-correlation function 
C{t), and then remains approximately constant for an 
interval at intermediate time scales, i.e. a "plateau" is 
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formed at those times. Finally, a renewed increase takes 
place on long time scales which is completely governed by 
the Drude weight, or by the zero- frequency distribution of 
C{uj). Note that the plateau can not be seen, if zero- and 
finite-frequency parts are treated separately from each 
other. It is a feature which arises from a combination of 
both contributions. 

The above plateau remarkably is already visible for 
= 8. Moreover, its "height" does not change with N, 
while its "width" seems to increase gradually, see Fig. |4] 
(middle inset). The latter increase particularly appears 
to be plausible, because the Drude weight is commonly 
expected to vanish in the thermodynamic limit N ^ oo, 
as already outlined above. We hence make the educated 
guess that the plateau of I?(i) will be continued to, say, 
arbitrary long times, when only TV becomes sufficiently 
large, e.g. N — > oo. Then the height P/J w 0.60 directly 
determines the concrete value of the diffusion constant, 
of course. And indeed, the latter value for T> also is in ex- 
cellent agreement with the conductivity I?bath/>/ — 0.58 
in Ref. Il3^. as found therein from a non-equilibrium 
bath scenario. Although a respective figure is not shown 
for A = 1.6, such an agreement is additionally obtained 
with the results in Ref. [l^ namely V/J « 0.55 and 
^'bath/^ — 0.585 ± 0.020. (In these works I?bath is de- 
noted by K.) Those quantitative agreements support the 
correctness of our guess and the emergence of diffusive 
transport appears to be verified for those A which are 
close to 1.5. 

Independent from the above arguments regarding the 
long-time extrapolation for — > cx), we should mention 
that for iV = 20 the 2?(t)-curves have already converged 
for those times which are shorter than i « 8/ J, see Fig. 0] 
(middle inset). Due to Eq. ((T2)) . this time corresponds to 
a length scale W on the order of about 3 sites. Thus, at 
least for Vl^ < 3 the dynamics should be known. 

For larger anisotropy parameters A we find a similar 
situation: Again there is a clearly developing plateau at 
intermediate time scales, see Fig. S] (bottom) for the case 
A = 2.0. For the accessible system sizes, however, the 
final plateau height V does not seem to be reached yet, 
i.e., V can not be read off with the same accuracy, as done 
before for A « 1.5. For comparison we hence indicate in 
Fig. m (bottom inset) the quantity a^c/iP^'^J) — 0.56 
with the dc-conductivity (Jdc/iPJ) — 0.14 according to 
Ref. [l^, obtained therein from an analysis on the basis of 
the standard Green-Kubo formula. The good agreement 
with this result is noticeable. 



C. Anisotropy A = 1 

As already mentioned above, it is still controversial, 
whether or not the Drude weight is finite in the limit 
N ^ oo, when A becomes exactly 1.0. However, let us 
for the moment assume that the Drude weight is indeed 
finite in that limit. Then at the infinite time scale the 
diffusion constant I?(i) is completely governed by this 



non-zero Drude weight and increases linearly. At finite 
time scales, however, 'D{t) may still appear to be almost 
constant and feature a plateau like the one for the case 
A = 1.5. In principle such a plateau can be very wide, 
if only its height is large in comparison with the Drude 
weight D such that the finiteness of D is by itself not 
crucial in this context, cf. Sec. HIl 

We hence show in Fig. 0] (top) the diffusion constant 
I?(t) for the special case A = 1. The initial short-time 
increase directly passes into the final long-time increase, 
i.e in between a horizontal line is not observable for the 
accessible system sizes. Nevertheless, we can not exclude 
the possibility that a plateau with a height on the order of 
V/J « 1.0 will eventually develop, when the system size 
is further increased, cf. Fig. [3] (top inset). We may thus 
compare this value for D with the extrapolated Drude 
weight D/J2 w 0.025 in Ref. for example. (The 
latter extrapolated D is about 1/4 of the long-time slope 
for N = 20, cf. Figs. [Hand [51) By the use of the rough 
estimation D\t2 — til/V ^ 1 we directly obtain that 
an interval with X'(i) « T) can not be wider than some 
time scale |t2 — ti\ <C 4:0/ J, respectively length scale 
\W2 — Wi\ -C 9. Thus, if the above extrapolation for 
the Drude weight was indeed correct, diffusion would be 
restricted to very few sites only. Contrary, if the Drude 
weight was zero (or a finite, irrelevantly small number), 
diffusion could occur for a possible arbitrary number of 
sites, of course. 



V. HEISENBERG CHAIN WITHIN AN 
ALTERNATING FIELD 

In this Section we are going to concentrate exclusively 
on the pure Heisenberg chain, i.e. A = 1, and investi- 
gate another modification of the Hamiltonian ([^T|) which 
may lead to diffusive transport. Such a modification cer- 
tainly is the incorporation of disorder, i.e. coupling con- 
stants or field strengths which vary randomly from one 
site to another.2iiSi2^iSii^. But the investigation of dis- 
order is numerically rather challenging and conceptually 
more subtle due to the localization phenomenon, even for 
a non-interacting single-particle syste m^^i^^ . Instead we 
will consider another scenario without disorder, where 
the Heisenberg chain is exposed to a strictly alternat- 
ing field. This scenario is concretely described by the 
Hamiltonian pTjl and an additional Zeeman term of the 
forniHJ^ 

H.^t^±S^S-. ,23, 

where 6B is the deviation from the mean B. As already 
done before, we may set B = for simplicity. 

The presence of the Zeeman term (|23|) does not affect 
the commutation of H and S^, i.e. the Hamiltonian can 
be diagonalized within decoupled M-subspaces, too. As 
long as N is even, there also is translational invariance 
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FIG. 6: The diffusion constant ©(t), as given by Eq. p3|) . for 
magnetization transport (spin transport) in the Ifeisenberg 
chain within an alternating field in the high temperature limit 
(T = cxa). Insets zoom in D(i) at short t. The curves are 
evaluated numerically (by the use of exact diagonalization) 
for chain lengths A'' = 10, 12, . . ., 18 (arrows) and for field 
strengths SB /J = 0.5 (top) as well as 0.75 (bottom). The 
crosses (insets) indicate the conductivities ©bath/ J = 0.323 ± 
0.010 and 0.15 ± 0.08 in Ref. [l^ . as obtained from a non- 
equilibrium bath scenario. 



(w.r.t two sites) such that the problem can be reduced 
further by a factor N/2. (There is no mirror symmetry 
for even N .) 

Since still represents a strictly conserved quantity, 
the present model allows to investigate the transport of 
magnetization w.r.t. the field parameter 5B. Note that 
the respective current and factor are identical to 
those in Sec. IIVI When SB is increased from zero and 
becomes comparative with the coupling strength J, the 
model undergoes a transition to quantum chaos^^,, i.e. one 
may assume that the latter transition already is a first 
pointer towards the onset of diffusion. Even though it 
is entirely independent from those assumptions, we start 
our investigation with 5B / J = 0.5. 

For 6B / J = 0.5 the curve for the diffusion constant 
'D{t) has changed from Fig. [3] (top) into Fig. [HI (top). In 
particular the initial short-time increase does not directly 
pass into the final long-time increase any more. Instead 
there is an oscillation in between such that, at the first 
view, the situation appears to be much more complicated 
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FIG. 7: The long-time increase D of the diffusion constant 
©(t), the Drude weight, for magnetization transport (spin 
transport) in the Heisenberg chain within an alternating field 
at high temperatures (T = oo). The symbols are evaluated 
numerically (by the use of exact diagonalization) for chain 
lengths iV < 18 and for field strengths SB/J < 0.75. 



than a simple horizontal line. However, for TV = 18 the 
2?(i)-curve has already converged until the end of this 
oscillation, see Fig. [5] (top inset). Exactly at this position 
the curve seems to gradually develop a plateau with a 
height on the order of P « 0.42 J. The plateau becomes 
visible, since the Drude weight rapidly decreases with 
iV, i.e. it is about one order of magnitude smaller for 
TV = 18 than for iV = 8, cf. Fig. [7] (circles). Note that the 
Drude weight does not fulfill a simple 1/iV-dependence 
and decays faster than a power-law. Thus, there is no 
need to suppose a finite and relevant Drude weight in 
the limit N ^ oo. 

Remarkably, the above suggested V for 6B/J — 0.5 
agrees well with the conductivity Pbath/ A = 0.323±0.010 
from Ref. [l6l . as found therein from a non-equilibrium 
bath scenario. The deviation from this conductivity is 
on the order of 30% solely. However, the small deviation 
may be explained as follows: In Ref. [l^ the conductiv- 
ity is evaluated by the use of the Lindblad equation, i.e., 
it is extracted from the steady state of a finite chain 
which is at both ends weakly coupled to baths at differ- 
ent temperatures. This weak bath-coupling is essential 
for the validity of the Lindblad approach. But, for a finite 
chain, a too weak bath-coupling may also yield a value 
2?bath which is lower than the correct value, say, V, see 
Ref. fsil. Since in Ref. [TBI the independence of the conduc- 
tivity from the bath-coupling strength is not analyzed in 
detail, we suppose that the above deviation results from 
a slightly too weak coupling of chain and reservoirs. 

For larger field parameters 6B suggestions are less re- 
liable, since the initial oscillations of the T>{t)-cmves be- 
come much more pronounced, as already visible in Fig. [6] 
(bottom). Even though there are oscillations, T>{t) still 
remains strictly positive, i.e. those oscillations probably 
are no pointer towards insulating behavior. 
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Note that definite conclusions can not be made for the 
parameter regime SB/ J < 0.5, since the curves for T>{t) 
continuously change from Fig. [4] (top) into Fig. [6] (top), 
i.e. one basically is concerned with the problem that the 
Drude weight governs the overall shape of the curve, at 
least for accessible system sizes. 

VI. SUMMARY AND CONCLUSION 

In the present paper we have investigated transport 
in several translationally invariant spin- 1/2 chains in the 
special limit of high temperatures. We have concretely 
considered spin transport in the anisotropic Heisenberg 
chain, in the pure Heisenberg chain within an alternating 
field, and energy transport in an Ising chain which is 
exposed to a tilted field. 

To this end we have firstly reviewed on a recently de- 
rived connection between the evolution of the variance 
of some "typical" inhomogcneous non-equilibrium den- 
sity and the current auto-correlation function at finite 
times'^^. In the limit of infinitely long times this connec- 
tion was shown to yield a generalized Einstein relation 
which relates the diffusion constant (in the absence of 
any external force) to the dc-conductivity (as the lin- 
ear response coefficient in the presence of an external 
force, i.e. as evaluated from the standard Green-Kubo 
formulcfc^i^S) . However, we have additionally demon- 
strated that the great advantage of the above connection 
is given by its direct applicability at finite times and for 
finite systems, e.g. at short times interesting signatures 
of an infinitely large system may be extractable for a 
system with an accessible size. 

By means of numerically exact diagonalization we have 
indeed observed strong indications for diffusive behavior 
in the considered spin chains for a range of model pa- 
rameters. Moreover, the suggested diffusion constants 
have been found to be in quantitative agreement with 
recent results on diffusion coefficients which were ob- 
tained for the same spin chains from numerically in- 
volved investigations of non-equilibrium bath scenarios 
in Refs. [Ulliillii. 



Amongst all those and our findings at high tempera- 
tures the emergence of diffusive transport of magnetiza- 
tion appears to be verified in the anisotropic Heisenberg 
chain with the anisotropy parameter A = l.S ^'^i^^ , de- 
spite the integrability of the model. It is known that 
the onset of quantum chaos is not a sufficient condition 
for diffusion^' but the latter result suggests also that 
non-integrability is not necessary at all. 

Unfortunately, it is still an open question whether spin 
transport in the pure Heisenberg chain is diffusive or bal- 
listic. Simply by the use of numerically exact diagonal- 
ization we were not able to reach those system sizes which 
are required for any definite conclusion on that question. 

Therefore the next step certainly is the application 
of approximate methods in order to obtain the current 
auto-correlation function at finite times, either numeri- 
cal ones, e.g. Suzuki- Trotter decompositions^!^, or ana- 
lytical ones, e.g. projection operator technique o^^i'^^'^^i^^ 
or moment method o^^i^^ . Approximative methods will 
also be in indispensable, when the investigation is ex- 
tended from one-dimensional spin chains to, e.g. more- 
dimensional spin lattices. 

Apart from the above methodic details, a physically 
interesting question is the dependence of transport on 
temperature, of course. Since the connection between the 
evolution of the variance and the current auto-correlation 
function is not restricted to infinite temperature, further 
investigations may be done in this direction, too. How- 
ever, we expect that for low temperatures the conver- 
gence of the diffusion constant at finite times is much 
slower, when the system size is increased. Thus, for low 
temperatures conclusions on the basis of a finite system 
may be less reliable. 
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